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q^. In this paper we investigate how the phase diagram of a U(l) symmetric Higgs- 

Yukawa system depends on the scalar self coupling A. The phase diagram of simi- 
lar models with continuous symmetry were extensively studied in the infinite scalar 
self coupling A = oo limit. Recent analytical and numerical calculations at zero self 
coupling showed qualitatively different phase diagram, raising the question of the A 
dependence of the phase diagram. Here we use analytical (large Nf, perturbative and 
mean field) approximations as well as numerical simulations to investigate the system. 
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1 Introduction 



The non-perturbative studies of the Higgs- Yukawa systems were motivated by the ever in- 
creasing top quark mass limit and the triviality problem of the Standard Model. The phase 
diagram investigations of four dimensional Higgs- Yukawa models produced several surprising 
non-perturbative results in recent years JT|. 

In the weak Yukawa coupling region a Higgs- Yukawa model with continuous symmetry 
can be in a ferromagnetic (FM) phase, a symmetric (SYM) phase or an antiferromagnetic 
(AFM) phase, depending on the scalar hopping parameter (k) values. At large positive 
k values the continuous symmetry is spontaneously broken and the system is in the FM 
phase. The model has a massive scalar particle, one or more massless Goldstone bosons and 
fermions with mass generated via the spontaneous symmetry breaking of the scalar field. 
For small k values the system is in the SYM phase, in which the theory contains degenerate 
massive scalars and massless fermions (assuming there are no bare fermion mass terms in 
the action). The FM and SYM phases are separated by a second order phase transition line. 
The critical behavior of the model along this FM-SYM phase transition line is expected to 
be governed by the perturbative Gaussian fixed point, where both the scalar and Yukawa 
couplings are marginally irrelevant. In the infinite cut-off limit the renormalized couplings 
vanish, y R = 0, X R = 0. As the pure scalar model is "trivial", i.e. it has no other fixed point 
than the perturbative one, similar behavior is expected even for strong scalar coupling as 
long as the Yukawa interaction is weak. At large negative k values the model is in the AFM 
phase. It is generally believed that this phase is separated from the symmetric phase by a 
second order phase transition line the same way as the ferromagnetic and symmetric phases 
are in the positive hopping parameter region. 

For large Yukawa couplings the situation is quite different. Quenched and unquenched 
Monte Carlo simulations of the model with naive fermions, supported by strong Yukawa 
coupling expansion, revealed the existence of non-perturbative symmetric and broken phases 
in the large Yukawa coupling region. These phases are separated by a second order phase 
transition line where the fermions of the model have large (at the order of the cut-off) mass, 
the fermions decouple in the continuum limit leaving a non-interacting scalar theory behind. 

Several numerical calculations investigated Higgs- Yukawa systems recently JT|. Apart 
from studies of the discrete symmetry model 0, all considered the limit of infinite scalar 
coupling. The simulations for both the U(l) and SU(2) symmetric systems observed the 
phase diagrams in agreement with the above descriptions. With naive fermions the phase di- 
agrams showed perturbative and non-perturbative (strong Yukawa coupling) SYM, FM and 
AFM phases. In addition a ferrimagnetic (FI) phase was found numerically at the interme- 
diate Yukawa coupling values, where both the magnetization and staggered magnetization 
are finite ||. All phase boundaries were claimed to be second order. There exists a special 
point in the phase diagram where three phases, symmetric, ferromagnetic and ferrimagnetic 
coexist. It was speculated that this point could be a non-trivial fixed point where the criti- 
cal behavior of the system might change. This scenario however could not be confirmed by 
Monte Carlo simulations |4j . It was found that the critical behavior of the system very close 
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to this point is still consistent with the perturbative predictions. 



The vanishing scalar coupling A limit received attention last year when it was shown in 
the large fermion number (Nf) limit that the Higgs- Yukawa model is equivalent (up to in- 
verse cut-off corrections) to the generalized Nambu-Jona-Lasinio type models [[J. The phase 
diagram was calculated exactly in the large Nf limit for arbitrary Yukawa and vanishing 
scalar couplings. Numerical simulations agreed surprisingly well with the large Nf predic- 
tions even for Nf = 2. However this phase diagram turned out to be very different from the 
infinite scalar coupling case. A strong first order phase transition line was observed at A = 0. 
It was found that the FM to SYM phase transition line, which is relevant to the Standard 
Model physics, stops at the first oder phase transition line. There is no sign for non-trivial 
critical behavior along this line. All the points on it belong to the attractive domain of the 
Gaussian fixed point. The FI phase, observed at A = oo ||, does not exist at A = 0. 

This paper is our first attempt to understand how the phase diagram and the critical 
properties of the U(l) invariant Higgs- Yukawa model changes from A = to the A = oo limit. 
In section 2 we study the Higgs- Yukawa model under various theoretical approximations and 
a typical phase diagram for small A is discussed. In section 3 we present our numerical results 
for A < 1 and compare them with the theoretical expectations. In section 4 we conclude this 
study and discuss the A = oo limit briefly. 



2 Analytical Considerations 
2.1 The model 

The lattice action for the U(l) chiral invariant Higgs- Yukawa model with naive fermions is 
defined as 

S = S f + S H . (1) 
The fermion part of the action Sf is given by 

S f = J2Mx)M(x,z)^(z) , i = l,2,...,N f /2 . (2) 

In eqn. (^|) the fermion matrix may be written as 

M(x, z) = J2lv [<W, 2 - tf*-A»,J + V [0i (x) + ijsfaix)) S x>z , (3) 

where 7 M , 75 are the Hermitian Dirac matrices and y stands for the Yukawa coupling. The 
scalar part of the action Sh in eqn. (|l|) is given by 



S H = -K^2<f) a (x) [(f> a (x +fi) + (j) a {x -y)\ + J2^ 2 a { 



X) 



X,fl 



+ E A [0aW-l] 2 , a = 1,2 . (4) 
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In the numerical simulations, we use detM^M for the fermion matrix to keep the partition 
function positive definite. This is equivalent to including an extra fermion species in the 
action eqn. ([!]). 



2.2 Large Nf limit 

The model defined in eqns.(|T|-^) can be solved in the large fermion number Nf limit. The 
A = case was discussed in Ref ||, and now we include the quartic term in the analysis. 

It is convenient to consider a modified form of the scalar action 



+ X N (y 2 a (x) - N f ) 2 . 
The usual lattice action eqns. (|l|-§) is obtained by identifying 

k n = C 2 k , \ N = C 4 X , y N = Cy , 
where the factor C satisfies the equation 

C 4 - (1 - 2X N N f )C 2 - 2X N = 
The scalar field ip a (x) is related to the original field (j> a (x) by 

if a (x) = (f>a(x)/C . 



(5) 



(6) 



(7) 



(8) 

In the 1/Nf expansion the couplings y^ = ^J~Nf~y N , X N = Nf\ N are kept fixed to be 0(1). 
As Nf — > oo the relations in eqn. (|6]) simplify, giving 



K N 



1 - 2A 



, XN f 



X 



N 



N 



1 - 2Aat) 2 



Vn 



1 -2A 



(9) 



N 



In the large Nf limit the constant mode of the scalar field dominates the path integral which 
suggests the Ansatz 



i Pl (x) = JN f (a + (-l)^^b) ,<p 2 (x) = 



(10) 



where yiV/a and JNfb correspond to the magnetization and staggered magnetization, re- 
spectively. The effective potential at leading order is 



N 



]K N {a 2 - b 2 ) + (a 2 + b 2 ) + X N ((a 2 + b 2 - I) 2 + Aa 2 b 2 ) (11) 



1 1 wv los 



sin 2 k^ + y 2 N (a 2 - b 2 



Minimizing V e ff (a, 6) with respect to a, 6 can give in general four types of solutions: 



(1) Symmetric (SYM) solution: a = 0, 6 = 0. 

(2) Ferromagnetic (FM) solution: a ^ 0, 6 = 0. The magnetization a is determined by 
the equation of state 



1 dV f 



eff 



N f da 

where the integral / is defined by 



= 2a(l - 8k n + 2(a 2 - 1)X N - 2y 2 N I(y 2 N a 2 )) = , (12) 

b=0 



r d 4 k 1 

I[X) ~ J {27xYY^sm 2 k, + x • (13) 

(3)Antiferromagnetic (AFM) solution: a = 0, 6^0. The staggered magnetization 6 is 
given by 

1 dV eff 



N f db 



26(1 + 8k n + 2(b 2 - 1)Ajv + 2y 2 N I(-y 2 N b 2 )) = . (14) 



a=0 



(4) Ferrimagnetic (FI) solution: a ^ 0,& ^ 0. It is straightforward to show that this 
solution can exist only when A^ > 1/2. However, the relation given in equation (0) becomes 
invalid for \ N > 1/2. Actually it is easy to find out from equation (|7|) that A at > 1/2 
corresponds to the intermediate and strong A region for the original lattice action eqn. ([I]), 
which we will investigate in a future study. Thus the large Nf calculation indicates that at 
least there is no FI phase in the weak A region. 

In some parameter range several solutions may coexist. It is a simple numerical exercise 
to evaluate the effective potential and find the solution that gives the absolute minimum. 

Although the above approach is exact in the limit Nf — > oo with fixed Aat and yN, it 
only describes the small A(= 0{1/Nfj) and y(= 0(l/-^iV/)) region of the original model 
eqn. ([I]). For large y values another type of large Nf expansion is possible 0. We start with 
the action eqns. ([!]) and (H), only now assuming that y N ~ 0(1) and k n ~ 0(1/Nf). After 
integrating out the fermions, we get an effective action for the scalar variables 

S eff = S H -^tr log MMt , (15) 

where Sh and M are given by eqns. (|5|) and (0), respectively. One may expand the fermion 
determinant in powers of 1/Nf. The leading term in 1/Nf fixes the amplitude of (f a {x) 



<p a {x) = Vo° a (x) , 1 + 2A^(^ - 1) - = , (16) 

Nf <p 2 Q 

where a a (x) is a two-component field with unit length. At next to leading order in 1/Nf, the 



effective action eqn. (|T^) becomes, up to an additive constant, an effective 4-dimensional 
XY-model 

S ef f = -K ef f ^2 a a (x) [a a (x + n) + a a (x - //)] . (17) 

X,fl 
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The effective hopping parameter K e ff is given by 



K ^ = KN "> 1 + 2yM ■ (18) 

The XY-model is known to have second order phase transitions at K e ff ~ 0.15 (FM-SYM) 
and n e ff ~ —0.15 (AFM-SYM). Thus eqn. fll%|) predicts the existence of two second order 
phase transition lines for large y^. 

Combining the above two large Nf expansion results, we plot the phase diagram for fixed 
Ajv(= 0.1) in Fig. 1. We find the following: 



(1) An FM-SYM second order phase transition line (AB). It is described by eqn. fll2|) in 
the limit a —>■ 0. 

(2) An AFM-SYM second order phase transition line (CD). It is given by eqn. (|14]) in 
the limit b — > 0. 

(3) A first order phase transition line (EDBF). On this line the effective potential has 
an AFM minimum in coexistence with either another AFM, SYM, or FM minimum. The 
staggered magnetization is discontinuous along the line DE but it is finite on both sides of 
the phase transition line. The discontinuity decreases from D to E and becomes zero at 
the point E. Therefore the AFM phases in the small and large y^ regions are analytically 
connected. For increasing \^ values the position of the DE segment moves to larger y^ and 
the end point E moves to more negative kn value. As Xn — > the line CD disappears and 
the point E coincides with the point C as given in Ref B. 

(4) A second order FM-SYM phase transition line (GH) given by eqn. setting 
K ef f = 0.15. 

(5) A second order AFM-SYM phase transition line (IJ) given by eqn. (|D|) setting 

K eff = -O- 15 - 

The region between F,K,I,G is outside the validity of both large Nf expansions. The 
order of phase transitions for line FK,KI, and KG can be only determined by numerical 
simulations. The second order FM-SYM transition line AB ends on a first order phase 
transition line but the first order line does not get critical at this point. The whole AB line, 
including the endpoint B, is in the domain of attraction of the Gaussian fixed point. No new 
non-trivial fixed point is found in this region. 

Although the above phase diagram is obtained in the large Nf expansions, the bare 
perturbation calculation in the small y region predicts the same structure for finite Nf as 
shown in Fig. 2a. 
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2.3 Mean field calculation 



When A is relatively large, both large Nj expansions and the bare perturbation calculation 
become invalid. For a qualitative phase diagram, one may use mean field calculations. Here 
we follow the well know |7|, j|, |6], || saddle point type mean field approximation. 



Because of the U(l) chiral symmetry, one may choose the Ansatz for the saddle point 

(19) 



fa(x) = a + (-l)^ Xfl b, 2 (x) = O, 
hi(x) = h + (-l)^ x »h st , h 2 (x)=0 



where hi(x),h,2(x) are the auxiliary fields introduced in the mean field calculation. The 
saddle point conditions define h and h s t as implicit functions of the magnetization a and the 
staggered magnetization 6 



a + -u'(h + h st ) + ~u'(h - h st ) = , 

1 l 

b + -u'(h + h st ) - -u'(h - hst) = . 



where u(x) is a function defined by 



exp{-u(x)} = J dpdOexp {—V(p) — px cos 6} 



and V(p) =p 2 + A(p 2 - l) 2 . 

With this Ansatz, the effective potential of the system has the form 



- V MF {a, b) = 8K{a 2 - b 2 ) + 2N f J j^yj™ 



sin 2 k^ + y 2 (a 2 — b 2 



+ ha + h st b + ^u(h + h st ) + ^u(h - h st ) . 

Minimizing Vmf{o-, b) with respect to a, b gives three possible phases 

(1) SYM phase: a = 0, h = 0, 6 = 0, h st = 0. 

(2) FM phase: a ^ 0, h ^ 0, 6 = 0, h st — 0. a and h satisfy the equations 

a + u'{h) = , 
d A k 1 



16m + AN f y 2 a J 



h = 



(2tt) 4 sin 2 ^ + y 2 a 2 

(3) AFM phase: a = 0, h — 0, 6 ^ 0, h st ^ 0. 6 and /i st satisfy the equations 

6 + u'(h st ) = , 



(20) 
(21) 

(22) 



(23) 



(24) 
(25) 



(26) 
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-16.6-4^/^ — l—^ + k, t = 0. (27) 

Using the fact that u'(x) is an odd and monotonically increasing function of x, one can 
show that the FI (a 7^ 0, b 7^ 0) solution is excluded in the mean field approximation. 

Similar to the large Nf calculation, in some parameter range several solutions may coexist. 
Again one needs to evaluate the effective potential and find the solution that gives the 
absolute minimum. 

The above mean field calculation breaks down in the large y region due to the following 
simple reasoning: When y is large, the fermion determinant can be expanded and to leading 
order in 1/y 2 we have an effective action 



S[4>] = -k E [p( x )p( x + lA a a{x)a a {x + p) + {p -> -p)} (28) 

X,fJ, 

+ E(P^) 2 : + ^lp(x) 2 "if -2N f \ogp 2 (x)} , 

x 

where the radial and angular notation for the field is used 

0i (sc) = p(x)ai(x), a^x) = cos(9(x)) , (29) 
02(e) = p(x)a 2 (x), a 2 (x) = sin(6(x)) . 

If we stay on the k = axis, the system becomes a collection of uncorrelated rotators and 
the U(l) chiral symmetry will be unbroken. An expansion around the k = axis will have 
a finite radius of convergence. Thus by analytical continuation we expect the system to be 
in the symmetric phase in the large y region around k = 0. However, this symmetric phase 
in the large y region is not predicted by the mean field calculation. 

We comment here that for the Higgs- Yukawa model the saddle point type mean field 
approximation is not equivalent to the variational type mean field approximation. With 
the saddle point type approximation, one is not guaranteed to get an upper bound on the 
free energy. If the fluctuation around the saddle point is large, the true free energy may be 
completely different from the saddle point estimate. In contrast, the variational type mean 
field approximation gives the rigorous upper bound of the free energy. Unfortunately, the 
variational calculation can not be completed without further approximation (small or large 
y expansion 0) for the Higgs- Yukawa model. 

A different type of mean field approximation can be performed in the large y region. If 
we expand the fermion determinant to next leading order in 1/y 2 , the action becomes 



s = -E 



X,fl 



N f 

Kp(x)p(x + p) + 



2y 2 p(x)p(x + p) 
+ E{p(^) 2 + ^PW 2 -l] 2 -2iV / logp 2 (x)} . 



(T a (x)a a (x + p) + (p -> -p) \ (30) 
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The radial interaction part may be approximated by its mean field value 

c n ^ - jo aft cm 

Jo°° dppe-^ 2 + x (p 2 -^ 2 - m f lo sp 2 ] ' { ' 

and the action eqn. (30) becomes an effective action for the XY-model 

S ef f = -K ef f ^2 a a (x) [a a (x + p) + a a (x - p)] , (32) 

x,n 

where the effective hopping parameter n e ff is given by 

, Nf 1 , 

K eff = K< P > 2 < - > 2 • (33) 

This approximation can be justified as a leading order 1/y 2 expansion with the assumption 
that k ~ 0(l/y 2 ). The XY-model has second order phase transitions at n e ff ~ 0.15 (FM- 
SYM) and n e ff ~ —0.15 (AFM-SYM). Thus for the original Higgs- Yukawa model we get the 
second order phase transition lines at 

k < p > 2 +^4 < - > 2 ^ ±0.15 . (34) 
V P 



The results of this mean field calculation are plotted in Fig. 2a and 2b in the large y 
region. We want to mention that for small A and y the mean field results are compatible 
with the bare perturbation calculations which are plotted in Fig. 2a. 

At A = 1 both large Nf approximation and perturbation theory break down. Fig. 2b 
shows the phase diagram predicted by the weak and strong y mean field calculation at A = 1. 



3 Numerical results 



We have performed numerical simulations at A = 0.0156, Nf = 2, 10 and A = 0.1, 1.0, Nf 



2. The Hybrid Monte Carlo method [1C] was used for the dynamical fermion simulations. 
Each molecular dynamics trajectory consists of 10 steps with step size chosen such that 
the acceptance rate is around 80%. To decide the order of the phase transition, we looked 
for hysteresis effects in the thermocycles. For each data point in the thermocycle about 20 
trajectories are used as warmup and 100-200 trajectories are used in the measurement. 

The magnetization v is defined as 

V=<\f¥a>, 4>a = T jJ2M X ) , (35) 



and the staggered magnetization v st is 



Vst =< 



r ?st,a >, 4>st,a = -^£(-l)^0 a (:r) , (36) 
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where L is the linear size of the lattice, v and v s t are measured and used as the order 
parameters to determine the phase. The measurements are done on 4 4 lattices. A 4 4 lattice 
is certainly not sufficient to distinguish second order and weakly first order transitions. 
However, the combination and the agreement of the numerical and analytical results allow 
us to determine the order of the phase transitions reliably. 

At A = 0.0156, we have simulation results for both Nf = 2 and 10. For Nf = 10, the 
data can be compared directly with the large Nf calculation. We find complete agreement 
with the theoretical predictions as shown in Fig. 1. Along the line of the first order phase 
transition where the hysteresis effects are very strong, it is not easy to determine exactly 
the position of the phase transition. But the amplitude of v and v st agree very well with 
the large Nf prediction. Thus the phase transition line predicted by the large Nf expansion 
should be reliable. For large y, the SYM region is quite narrow for Nf = 10, making it 
difficult to establish this region numerically. However, the existing results || at A = and 
our simulation results presented in the next paragraph at Nf = 2 should be sufficiently 
convincing that the SYM phase indeed exist in this region. 

The phase diagram for A = 0.0156, Nf = 2 is plotted in Fig. 2a. The transition points 
agree well with both the bare perturbation calculation and the large Nf expansion in the 
small y region. The agreement with the large Nf calculation is probably due to the fact that 
the effective fermion flavor number around second order phase transition lines is 32 because of 
the lattice doubling effect. Up to A = 1 the picture remains qualitatively the same. The phase 
diagram for A = 1, Nf = 2 is shown in Fig. 2b. As A = 1 is already a coupling of intermediate 
strength, neither bare perturbation theory nor the large Nf approximations agree with the 
data points. However, the mean field approximation -from which the curves are plotted in 
the figure- shows reasonable agreement with the simulation result. In particular, there is no 
indication for an FI phase or a non-trivial fixed point at any of these coupling values. 

In Fig. 3 we give examples for thermocycles along the first order phase transition line 
EDBF in Fig. 2 (see Fig. 1 for notations) for A = 0.0156, Nf = 2. The data is compared 
to the bare perturbation theory (dotted lines). Fig. 3a corresponds to an AFM-AFM [] 
transition, Fig. 3b to a SYM-AFM and Fig. 3c to an FM-AFM transition. In all cases 
v st is plotted. All thermocycles show hysteresis effects in agreement with our theoretical 
expectations. There is no sign of the first order line becoming critical at the end point B of 
the FM-SYM second order phase transition line. 



4 Conclusion and discussions 

In this paper we have investigated the phase diagram of a U(1)<8>U(1) Higgs- Yukawa model 
with fluctuating length Higgs field. The results obtained by various analytical methods 
and numerical simulations show, up to the Higgs self coupling A = 1, a phase diagram as 

^^Note that this figure also contains the continuous phase transition from the symmetric to the antiferro- 
magnetic phase at k « —0.16 
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summarized in Fig. 1. This picture is quite different from the published results at A = oo. 
We observe strong first order phase transitions from the antiferromagnetic phase to the 
ferromagnetic, symmetric and antiferromagnetic phases. The second order line between the 
perturbative ferromagnetic and the symmetric phases, that can be relevant for the Standard 
Model, ends on this first order phase transition line. At this end point the first order line 
does not become critical, no new non-trivial critical behavior is expected, the whole second 
order line is in the domain of attraction of the Gaussian fixed point. In addition, we do not 
find a ferrimagnetic phase up to A = 1. 

In order to understand the critical behavior of the Higgs- Yukawa model one has to find 
out if the above picture for the phase diagram remains the same as one increases A further. 
If this is indeed the case, it would provide a natural explanation of why the critical indices 
do not change along the second order line at A = oo0. Another possibility would be that 
the phase diagram changes qualitatively at some strong A value and eventually merges into 
the published result at A = oo . The mean field picture given in section 2 remains 
unchanged up to A = oo. However, it is unclear how reliable these approximations really 
are. To determine the critical properties of the large and infinite A systems requires much 
more computer time, or different analytical approaches. This work is in progress and the 
result will be reported in a future publication. 
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Figure 1: The phase diagram in the large Nf limit at Xn = 0.1 and Nf = 10. The MC data 
are indicated by circles where the solid symbol denotes second and the open symbols first 
order phase transitions. The solid and dashed lines are the results from the 1/Nf expansions, 
where the solid lines represent second order and the dashed line first order phase transitions. 
In the middle of the phase diagram, where the 1/Nf expansions break down, the dotted lines 
indicate how the phase transition might continue. 



Figure 2: Phase diagram at A = 0.0156 (a) and A = 1.0 (b), both with Nf = 2. Here as 
in Fig.l solid symbols denote second and open symbols first order phase transitions. In the 
small y region the solid and dashed lines are obtained by (a) bare perturbation calculation; 
(b) mean field calculation. In both (a) and (b) the lines in the large y region are obtained 
from the mean field theory given in the second part of the section 2.3. Solid lines represent 
second order and dashed lines first order phase transitions. The dotted lines only indicate a 
possibility how the phase transition lines may continue. 



Figure 3: Hysteresis effects for the staggered magnetization v st are shown along the first 
order phase transition line EDBF (see Fig.l). The data are taken at A = 0.0156 and Nf = 2. 
The solid symbols represent the first half of the thermocycle and the open ones the way back. 
The solid lines are only connecting the data points to guide the eye. The dotted lines are the 
results from perturbation theory, (a) a point taken between D and E (AFM-AFM) (Note 
that in addition to the AFM-AFM phase transition also the SYM-AFM phase transition at 
k —0.16 is shown.) (b) a point between D and B (SYM-AFM). (c) a point taken between 
B and F (FM-AFM). The discontinuity of the hysteresis loop becomes smaller closer to point 
E indicating a weaker first order phase transition. 
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